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Slow oscillations of the in-plane magnetoresistance are observed in the rare-earth tritellurides 
and proposed as an effective tool to determine the parameters of electronic structure in various 
strongly anisotropic quasi-two-dimensional compounds. These oscillations do not originate from the 
small Fermi surface pockets, as revealed usually by the Shubnikov-de-Haas oscillations, but from 
the entanglement of close frequencies due to a finite interlayer transfer integral tz, which allows to 
estimate its value. For TbTes and GdTes we obtain the estimate tz « 1 meV. 
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The angular and magnetic quantum oscillations 
(MQO) of magnetoresistance (MR) is a powerful tool 
to study electronic properties of various quasi-two- 
dimensional (Q2D) layered metallic compounds, such as 
organic metals (see, e.g., Refs. [H-l^ for reviews), cuprate 
and iron-based high-temperature superconductors (see, 
e.g., iQ), heterostructures fisj etc. 

The Fermi surface (FS) of Q2D metals is a cylinder 
with weak warping At^/Ep <C 1, where tz is the inter¬ 
layer transfer integral and Ep = ^ is the in-plane Fermi 
energy. The MQO with such FS have two close funda¬ 
mental frequencies Eq ± AE. In a magnetic field B = Bz 
perpendicular to the conducting layers Eq/B — pi/huJc 
and AE/B = 2tz/tiu>c, where huc = heBz/m*c is the 
separation between the Landau levels (LL). 

The standard 3D theory of galvanomagnetic properties 
IMl is valid only at tz and in the lowest order 

in the parameter huc/tz- This theory predicts several 
peculiarities of MR in Q2D metals: the angular magne¬ 
toresistance oscillations {AMRO)[Ii-[2H and the beats of 
the amplitude of MQO [l7| . One can even extract the fine 
details of the FS, such as its in-plane anisotropy [ 2 ^ and 
its harmonic expansion, [illii] from the angular depen¬ 
dence of MQO frequencies and from AMRO. For isotropic 
in-plane electron dispersion, AMRO can be described by 
the renormalization of the interlayer transfer integral: [2^ 


tz = tz {0) = tz (0) Jo {kpcLianO ), 


( 1 ) 


where Jo {x) is the Bessel’s function, pp = hkp is the 
in-plane Fermi momentum, and 6 is the angle between 
magnetic field B and the normal to the conducting layers. 

At tz fvjjc new qualitative features appear both in 
the monotonic and oscillating parts of MR. For example, 
the strong monotonic growth of inte rlay er MR Rzz{Bz) 
was observed in various Q2D metals 26l435l| and recently 


theoretically explained [35l - l38j . The oscillating part of 
interlayer MR at p ^ tz ^ hujc acquires the slow oscilla¬ 
tions [il, [ 3 ^ and the phase shift of beats, li^ These 
two effects are missed by the standard 3D theory [iMi 
because they appear in the higher orders in huc/tz- 

The slow oscillations qualitatiyely originate from the 
product of the oscillations with two close frequencies 
Eq ± AE, which giyes the oscillations with frequency 
2AF. Conductiyity, being a non-linear function of the 
oscillating electronic density of states (DoS) and of the 
diffusion coefficient, has slow oscillations with frequency 
2 AE, while magnetization, being a linear functional of 
DoS, does not show slow oscillations 34, The slow 
oscillations haye many interesting and useful features as 
compared to the fast quantum oscillations. First, they 
suryive at much higher temperature than MQO. Second, 
they are not sensitiye to the long-range disorder, which 
has a strong action on fast MQO, similar to the influ¬ 
ence of finite temperature. Third, the slow oscillations 
allow to measure the interlayer transfer integral tz. These 
features make the slow oscillations to be a useful tool to 
study the electronic properties of Q2D metals. Until now, 
the slow oscillations where inyestigated only for the in¬ 
terlayer conductiyity azz {B), when the current and the 
magnetic field are both applied perpendicularly to the 2D 
layers, and only in the organic compounds. In this paper, 
we study the slow oscillations of transyerse intralayer MR 
Ryy{B) both theoretically and experimentally in the non- 
organic rare-earth tritelluride Q2D layered compounds. 

The rare-earth tritelluride compounds RTe^ {R =Y, 
La, Ce, Nd, Sm, Gd, Tb, Ho, Dy, Er, Tm) haye the same 
orthorhombic structure {Cmcm) in the normal state. 
These systems exhibit a c-axis incommensurate charge- 
density waye (GDW) at high temperature through the 
whole R series that was recently a subject of intense stud¬ 
ies (dll - li^ . For the heayiest rare-earth elements, a second 
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FIG. 1: (color online) dR{B)/dB dependencies at T = 4.2 K 
for TbTea (a) and GdTea (c) demonstrating clear rapid oscil¬ 
lations which appear in the field range B > 2 T. Insets show 
magnetoresistance R{B). (b) and (d) - dR/dB as a function 
of inverse magnetic field, B~^, at low field region demonstrat¬ 
ing slow oscillations in TbTes and GdTes correspondingly. 
Insets show the Fourier transform of slow oscillations. 


a-axis CDW occurs at low temperature. In addition to 
hosting incommensurate CDW, magnetic rare-earth ions 
exhibit closed-spaced magnetic phase transitions below 
10 K [13,13 leading to coexistence and competition of 
many ordered states at low temperature. Therefore, any 
information about the Fermi surface warping in the Q2D- 
compounds (determined by tz) is very important. For 
the possible observation of the slow oscillations the rare- 
earth tritellurides are very promising, because they have 
appropriate anisotropy and good metallic conductivity 
up to low temperatures. 

For experiments we choose GdTes and TbTes. Single 
crystals of these compounds were grown by a self-flux 
technique under purified argon atmosphere as described 
previously [3 ■ Thin single crystal samples with a thick¬ 
ness typically 0.1-0.3 ^m were prepared by micromechan¬ 
ical exfoliation of relatively thick crystals glued on a sap¬ 
phire substrate. The quality of selected crystals and 
the spatial arrangement of crystallographic axes were 
controlled by the X-ray diffraction. From high-quality 
[R{300K)/R{10K) > 100] untwinned single crystals we 
cut bridges with a width 50 — 80 ^m in well defined, 
namely [100] and [001], orientations. Contacts for elec¬ 
trical transport measurements in the four-probe config¬ 
uration have been prepared using gold evaporation and 
cold soldering by In. Magnetotransport measurements 
were performed at different orientations of the magnetic 
field in the field range up to 9 T using a superconducting 
solenoid. The field orientation was defined by the angle 9 
between the field direction and the normal to the highly 
conducting ac plane. 

In Figure [T] we plot the derivative of MR dR/dB as a 


function of magnetic field up to B = 8.2 T applied along 
6 -axis for TbTea (a) and GdTes (c) with the current ap¬ 
plied in the (a, c) plane at T = 4.2 K. For both measured 
compounds, at B > 2T the pronounced Shubnikov - de 
Haas (SdH) oscillations with a frequency F Ri 55 — 58 
T are observed in the derivative of MR, dR/dB. At 
high field {B > 7 T) new oscillations with higher fre¬ 
quency {F « 0.7 — 0.8 kT) appear, indicating the exis¬ 
tence of several types of pockets on the partially gapped 
FS, resulting from the CDW formation Insets in 

Fig. [T^,c show the corresponding MR curves indicating 
very weak amplitude of oscillations. In addition to rapid 
SdH oscillations, at low magnetic field (H < 2 T) the 
MR exhibits prominent slow oscillations with very low 
frequency Fgiow ^ 4 T. To demonstrate the periodicity 
of these oscillations in we plot the derivative dR/dB 

as a function of inverse magnetic field in Fig. [T]for TbTea 
(b) and for GdTes (d). Insets show the Fourier transform 
of dR/dB{l/B). 

In contrast to the usual SdH oscillations, the amplitude 
of which decreases rapidly as temperature increases, the 
slow oscillations of MR are observable up to T 40 K, as 
can be seen from Fig. where we show the temperature 
evolution of slow oscillations for GdTes and TbTes. This 
suggests that the observed slow oscillations originate not 
from small FS pockets, but from the FS warping due to 
tz, similarly to the slow oscillations of interlayer MR in 
the organic superconductor/3-(BEDT-TTF)2lBr2 34j. If 
so, the observed slow oscillations give an excellent oppor¬ 
tunity to find the value of tz and of the product kpd at 
low temperature in tritellurides TbTes and GdTes from 
the experimental data on the intralayer MR. 

According to Eq. (90.5) of Ref. 3 the intralayer 
conductivity at finite temperature is given by (3 




( 2 ) 


where the derivative of the Fermi distribution func¬ 
tion n'p{£) = —1/{4Tcosh^ [(e — ^)/2r]}, and the zero- 
temperature electron conductivity at energy e is 


Uyy^e) =e^g {e) Dy (e). 


(3) 


Here g (e) is the DoS and Dy (e) is the diffusion coeffi¬ 
cient of electrons along y-axis. It is convenient to use the 
harmonic expansion for the oscillating DoS g^e). Below 
we will need only the first terms in this harmonic series, 
which at finite tz ^ Jtuic are given by Sim HI 


9 (e) - 90 


/ 2'K£ 




( 4 ) 


where go = m*/irh^d is the DoS at the Fermi level in the 
absence of magnetic field per two spin components, Jo (cc) 
is the Bessel’s function, the Dingle factor 3 Rd ^ 
exp [—Trfc/wcTo]: A) = fi/2To is the electron mean free 
time without magnetic field, and Fo is the LL broadening. 
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FIG. 2: (color online) Temperature evolution of slow oscilla¬ 
tions for GdTea (a) and TbTea (b). 


where t (e) is the energy-dependent electron mean scat¬ 
tering time by impurities, and the angular brackets in 
Eq. (|ni) mean averaging over impurity scattering events. 
In the Born approximation, the mean scattering rate 

l/r (e) = 2TrniU^g (e), (7) 


where rii is the impurity concentration. This scattering 
rate has MQO, proportional to those of the DoS in Eq. 
dl]). The MQO of are, usually, weaker 

and in 3D metals they are neglected [d^. Then 

Dy{e) ^ Rl/2T{e) (X g{e). (8) 


However, in Q2D metals, when tz ~ fiwc, the MQO of 
can be of the same order as MQO of the DoS. 
Then instead of Eq. at Rjj <C 1 one has 


Dy {e) 


Do 


1 — 2a cos 





(9) 


where Dq sa R\/2tq, and the number a 1. Substitut¬ 
ing Eqs. ®,o and ® to Eq. @ after integration over 
e one obtains 

^ 1 + 2aJo" {^TTtz/hwc) Rl- (10) 

e^goDo 

_2(a + l)cos(g) 


The calculation of the diffusion coefficient Dy (e) is less 
trivial and requires to specify the model. At g ^ fnoc the 
quasi-classical approximation is applicable. In an ideal 
crystal in a magnetic field B the electrons move along 
the cyclotron orbits with a fixed center and the Larmor 
radius Rl = ppc/eBz- Without scattering the electron 
diffusion in the direction perpendicular to B is absent. 
Scattering by impurities changes the electronic states and 
leads to the electron diffusion. For simplicity, we consider 
only short-range impurities, described by the 5-function 
potential: Vt (r) = U5^ (r — ri). Scattering by impurities 
is elastic, i.e. it conserves the electron energy e, but the 
quantum numbers of electron states may change. The 
matrix element of impurity scattering is given by 

Tmm' = (d) U^m (u) , (5) 


where 'I'm (r) is the electron wave function in the state 
m. During each scattering, the typical change Ag = 
APxc/eBz of the mean electron coordinate yo perpen¬ 
dicular to B is of the order oi Rl, because for larger 
Ay ^ Rl the matrix element in Eq. ([^ is exponentially 
small because of small overlap of the electron wave func¬ 
tions T'm, {vi) 4'm in) ~ 4'm (d + ^v) (?'i) HI. The 

diffusion coefficient is approximately given by 


Dy (e) ~ {{Ayf"^ /2t (e) 


( 6 ) 


where the temperature damping factor of MQO is 

Rt = (2t:'^ k bT/ huj^ / smh{2-K‘^kBT/huj^ . (11) 


The slow oscillations, described by the first line of Eq. 
(nnj, are not damped by temperature within our model. 

Approximately, one can use the asymptotic expansion 
of the Bessel function in Eq. m for large values of 
the argument: Jo(x) ~ ^21'kx cos (cc — 7r/4) , x » 1. 
Then, after introducing the frequency of slow oscillations 
Fsiow = ^tzB/hwc, the first line in Eq. (fTOl) simplifies to 


vv 


\B) 


(F-goDo 


ahuJr 


'2ttF, 


slow 


\ B 


R 


D- 


( 12 ) 


In tilted magnetic field at constant \B\, ujc oc cos0 and 
tz changes according to Eq. O- Then the frequency of 
slow oscillations must depend on tilt angle 6 as: 


Fsiow (0) /Fsiow (0) = Jo {kpd tan 0) / cos {0). (13) 


To clarify this effect, we experimentally study the an¬ 
gular dependence of observed slow oscillations frequency. 
The evolution of the slow oscillations in GdTea with the 
change of the tilt angle 0 of magnetic field at T = 4.2 K 
is shown in Fig. [U where the derivative dR/dB is plotted 
as function of the perpendicular-to-layers component of 
magnetic field B± = Bcos{0). Note, that magnetic field 
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FIG. 3: (color online) Slow oscillations in GdTes at T = 4.2 
K for different tilt angles Q. Bx = Bcos{0). 



FIG. 4: (color online) (a) and (c): The frequency of slow os¬ 
cillations Fsiow as a function of tilt angle 9 at T = 4.2K for 
TbTea and GdTes. Solid curves show the function F{6) — 
F{0)/cos{9). (b) and (d): Angular dependence of the fre¬ 
quency Fsiowifi) of slow oscillations of intralayer magnetore- 
sistance in TbTes and GdTes correspondingly, multiplied by 
cos(0). The experimental data are shown by black filled 
squares, and the theoretical prediction according to Eq. (na) 
at kpd = Q.12 is shown by solid red lines. 


rotation in the {b-c) and (6-a) planes demonstrated the 
same results for TbTes. 

In Fig. 0] we show the 0-dependence of the frequency 
of slow oscillations at T = 4.2 K for TbTes (a) and for 
GdTes (c). The solid curves give the cosine dependence 
F{9) = F{Q)/ cos(0) typical for MQO. According to Eq. 
m, the angular dependence of the frequency Fsiow (0) of 
slow oscillations differs from this standard cosine depen¬ 
dence, especially at high tilt angle. In Fig. 0] (b) we plot 
the angular dependence of the product Fsiow{d) cos(0) in 
TbTes. If the origin of slow oscillations is due to small FS 


pockets, the product Fsiow{d) cos{6) would be indepen¬ 
dent of the tilt angle 0. The experimental data, shown 
by black filled squares, clearly indicate the deviation from 
the horizontal line. These experimental data can be rea¬ 
sonably fitted by Eq. (fT51) at kpd — 0.1, shown by solid 
red lines in Figs. 0] (b,d). 

From the angular dependence of the frequency 
Fsiow (0) of slow oscillations, shown in Fig. 01 we con¬ 
clude that these oscillations originate not from small FS 
pockets as usual SdH oscillations, but from the entan¬ 
glement of close frequencies due to a finite tz, similar to 
that observed in the interlayer MR in the organic metal 
/3-(BEDT-TTF)2lBr2 34|. The strong additional argu¬ 
ment in favor of this origin of the observed slow oscilla¬ 
tions is the very weak temperature dependence of their 
amplitude. To our knowledge, the data obtained are the 
first observation of such slow oscillations in the intralayer 
magnetotransport. 

The frequency of slow oscillations at 0 = 0 can be used 
to estimate the value of tz- According to Eq. (II2p . with 
the effective electron mass m* Ri O.lrue determined from 
the temperature dependence of the amplitude of SdH os¬ 
cillations 54|, and Fgiow ~ 3.5T (see Fig. 0] (b,d)), we 
obtain tz ~ ImeE. The small values of tz compared 
to the transfer integrals t|| ~ 2 eV along the chains 
and t± Ki 0.37 eV perpendicular to the chains, as ob¬ 
tained by the band structure calculations [42|, illustrate 
the quasi-2D character of these rare-earth tritellurides 
and justify that the dispersion along 6-axis is neglected 
in ARPES measurements. The value of tz is very im¬ 


portant for quantum corrections to conductivity [57 
and for other physical properties of strongly anisotropic 
compounds. 

The angular dependence of the frequency Fsiow (0) of 
slow oscillations allows also to determine the value of 
the Fermi momentum of the open FS pockets. Fitting 
experimental data on Fsiow (0) shown in Fig. 0] to Eq. 
(fO)l gives kpd ~ 0.11 for GdTes and kpd « 0.12 for 
TbTes. With d « 25A this gives kp Ri 4.5 • 10^cm~^. 
Such a small value of kp means that the FS reconstruc¬ 
tion due to GDW is very strong and leaves only very 
small ungapped FS pockets. These FS pockets are most 
probably elongated and directed differently, which gives 
different Fsiow (0) dependence. Their total contribution 
to the slow oscillations, being a sum of the contributions 
from the individual FS pockets, has a smeared angular 
dependence of the frequency of slow oscillations Fsiow (0) 
as compared to the case of only one elliptical FS pocket 
observed in /3-(BEDT-TTF)2lBr2 

To summarize, we report the hrst observation and 
qualitative theoretical description of the slow oscillations 
of intralayer magnetoresistance in quasi-2D metallic com¬ 
pounds. These slow oscillations allow to measure the in¬ 
terlayer transfer integral tz, which is hard to measure by 
any other ways. We obtain the value tz ^ ImeV in the 
rare-earth tritelluride compounds TbTes and GdTes. 
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